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Abstract 



o 



)rhis paper addresses the problem of blind separation of convolutive mixtures of BPSK and circular linearly modulated 
signals with unknown (and possibly different) baud rates and carrier frequencies. In previous works, we established 
|that the Constant Modulus Algorithm (CMA) is able to extract a source from a convolutive mixture of circular linearly 
modulated signals. We extend the analysis of the extraction capabilities of the CMA when the mixing also contains BPSK 
Signals. We prove that if the various source signals do not share any non zero cyclic frequency nor any non conjugate 
cyclic frequencies, the local minima of the constant modulus cost function are separating filters. Unfortunately, the 
minimization of the Godard cost function generally fails when considering BPSK signals that have the same rates and 
the same carrier frequencies. This failure is due to the existence of non-separating local minima of the Godard cost 
|function. In order to achieve the separation, we propose a simple modification of the Godard cost function which only 
requires knowledge of the BPSK sources frequency offsets at the receiver side. We provide various simulations of realistic 
^digital communications scenarios that support our theoretical statements. 
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1. Introduction 

] The blind source separation of convolutive mixtures of 
•linearly modulated signals has mainly been studied in the 
.case where the signals share the same known baud rate, 
and when the sampling frequency of the multivariate re- 
'ceived signal coincides with this baud-rate. In this context, 
to be referred to in the sequel as the stationary case, the 
|discrete-time received signal coincides with the output of 
an unknown MIMO filter driven by the sequences of sym- 
]bols sent by the various transmitters. In most cases, these 
sequences are independent and identically distributed, and 
^several methods have been proposed in order to extract 
'each of them from the observation (see e.g. Q, Q, 0, 
12 1, [l^) . The source separation problems that are en- 
countered in the context of passive listening are however 
more complicated because the transmitters are usually 
completely unknown to the receiver, and have no reason 
to transmit linearly modulated signals sharing the same 
baud-rates. It is therefore quite relevant to address the 
problem of blind separation of linearly modulated signals 
with unknown, and possibly different, baud rates. In this 
context, the received signal is sampled at any frequency 
satisfying the Shannon sampling theorem, so that the cor- 
responding discrete-time signal is cyclostationary with un- 
known cyclic frequencies. If the cyclic frequencies were 
known at the receiver side, it would be easy to generalize 
the usual blind source separation approaches based on the 



optimization of contrast functions depending on higher or- 
der cumulants. However, when the cyclic frequencies are 
unknown, it is impossible to consistently estimate the cu- 
mulants, a conceptual problem first remarked by Ferreol 
and Chevalier (Q) in the context of blind separation of in- 
stantaneous mixtures. An obvious approach would consist 
in estimating the unknown cyclic frequencies. However, 
this is a difficult task if the excess bandwidths of the trans- 
mitted signals arc low and if the duration of observation 
is not large enough. 

In contrast with the cumulants, the constant modulus 
cost function can be consistently estimated in the cyclosta- 
tionary context. In [loj . we considered only source signals 
that transmit second-order circular symbol sequences, and 
we have shown that in this case, to be referred to as the cir- 
cular case, the minimization of the Godard cost function 
allows to extract the sources using a defiation approach 
if their baud-rates are different one from another. If cer- 
tain baud rates coincide, sufficient conditions for the sep- 
aration have been established in (T^]. Although we have 
not been able to prove that separation is achieved in the 
most general case, all the simulations we have performed 
strongly suggest that the minimization of the Godard cost 
function is successful in the circular case. The purpose of 
this paper is to address this issue when in the non circular 
source signals, which will be referred to as the non circular 
case, and to show how the separation method based on the 
minimization of the CMA contrast function coupled with 
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a deflation approach can be adapted to this context. As 
in we only focus in this paper on the separation of the 
first source. 

In order to simplify the presentation of our results, we 
only consider the case where the non circular signals are 
BPSK signals. We begin by defining in section [5] the con- 
text of our study and giving a brief description of the con- 
sidered signals and criteria. In section |3] we prove that the 
Godard cost function is still successful if the sources do not 
share the same baud rates and the same carrier frequencies. 
We also prove, in section |H that contrary to the circular 
case, the minimization of the Godard cost function fails to 
separate 2 BPSK signals sharing the same baud rate and 
the same carrier frequency. We show that this is due to 
the existence of non separating local minima of the Godard 
cost function, toward which the minimization algorithms 
seem to converge quite often. Wc also show that it is possi- 
ble to modify the Godard cost function in order to achieve 
source separation of K non circular BPSK modulated sig- 
nals sharing the same known (or well estimated) carrier 
frequency. Section [S] briefly generalizes this result to more 
general mixtures. The new modified CMA algorithm needs 
the estimation of the carrier frequencies offsets of the non 
circular source signals, or equivalently the estimation of 
the "significant" non conjugate cyclic frequencies of the re- 
ceived signal. Fortunately, this is a much easier task than 
the estimation of baud rates, because the non conjugate 
cyclic correlation coefficients of the received signal at twice 
the frequency offsets are not affected by possible low ex- 
cess bandwidths of the source signals (see Numerical 
results are finally presented in section [51 

Notations: If (u„)„f:z is a discrete-time sequence, we de- 
note by < u„ > the time average operator defined as 



< u„ >- 



lim 

N^+cx, 2N + 



If X is a complex valued random variable, we denote by 
C4(x) its fourth order cumulant defined by cum{x,x* ,x,x*}. 
If (jc(n))„£z is a discrete-time cyclostationary sequence, we 
define, when it makes sense, the cyclo-correlation at cyclic- 
frequency a and time lag m: 



Vore 
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Vm € Z, R^^\m) =< E(x(n + OT)x(n)*e"2"™) > 



and the non conjugate cyclo-correlation at cyclic-frequency 
Qc and time lag m: 



' 2 2 



2iKna\ 



For an interval S, wc denote by 'F(S) the set of all func- 
tions fa{t) € 1?{W) such that 



Js 



,2mvl 



fa(y)dv 



In other words, a square integrable function is an cle- 
ment of "FiS) if and only if its Fourier transform /a(v) is 
zero outside S. 

2. Problem statement 

2.1. Assumptions 

We assume that K unknown transmitters send linearly 
modulated signals sharing the same frequency bandwidth. 
The receiver is equipped with a sensor of A^-arrays, and the 
corresponding A^-dimensional received signal is sampled at 
rate supposed to satisfy the Shannon sampling theorem. 
For any k, k — I, . . . ,K, the signal transmitted by source 
k is obtained by linearly modulating a unit variance zero 
mean i.i.d. sequence of symbols {akj,}„ez with a shaping 
filter ga,k 

Wc denote by Tk the symbol period of the source number 
k and wc consider a shaping filter of limited bandwidth 
[-^^, ^9^], where jk is the excess bandwidth factor, be- 
longing to [0, 1). The bandwidth of the complex envelope 
of transmitted signal k is then [--^^, 

In order to simplify the presentation of the results we 
make the following assumption: 

• the symbol sequence {flA-,„)„ez is either second order 
circular or corresponding to a BPSK constellation (i.e. 
equal to +1) for each k. 

The propagation channels between each transmitter and 
the receiver arc assumed to be frequency selective. More- 
over, the carrier frequencies of the various transmitted sig- 
nals of course do not coincide with the center frequency of 
the receive filter of the receiver. Hence, the contribution 
of each transmitted signal at the receiver side is corrupted 
by a frequency offset. The frequency offset associated to 
source k is denoted by Afk. 

We denote by yakiO the dimensional continuous-time 
signal representing the contribution of the transmitted sig- 
nal k to the received signal ya(t) which is to say, the signal 
that would be received if only transmitter k were active. 
Wc can then write ya,kit) as 



ya.k(t)^e'-'''^A'iK,k*Sa,k)(t) 



(1) 



For a wide-sense cyclostationary continuous-time ran- 
dom process {xa(t)),em wc denote by R'a,i'{T) and by R^a,c^(T) 
the cyclic correlation coefficient and respectively non con- 
jugate cyclic correlation coefficient at cyclic-frequency 
ttfl (respectively non conjugate cyclic frequency aa,c) and 
time lag t. 



where * represents the convolution operator and where ha^k 
is the dimensional channel impulse response between 
source k and the multiple-sensors receiver. The presence 
of the frequency offset shifts the bandwidth of the ya,k(0 
signal with a factor equal to Afk, thus making it coincide 
with the interval [-^-^ + Afk, ^ + Afk]. 



The continuous-time received signal (in the absence of 
noise) ya(t) = 2f=i ya,kiO is sampled at rate Tg which is 
supposed to verify 



— >max^ + |A/, 



(2) 



Under these assumptions, the A^-dimensional discrete-time 
received signal y{n) can be written as 

K r ^ ^ 

k=l V / / i=l 

(3) 

where for each k, stin) represents the sampled version of 
transmitted signal k, and where hjt(z) = Yjiez ^k.iz ' is the 
transfer function of the 1 -input / outputs discrete time 
equivalent channel between transmitter k and the receiver. 
Finally, dfk is defined as 6fi< = AfkTe- 

2.2. Expansion of the Godard cost function 

Due to the previously described context, each of the 
transmitted signals is cyclostationary and thus has a set 
of cyclic frequencies which are easily identified from the 
second order statistics of each signal. For all k, and for all 
T G K, the cyclic correlation function f — > E(ifl,jt(f +T)ia,i(f)*) 
and, for a BPSK signal, the non conjugate cyclic correla- 
tion function t — > E{sa,kit + T)sa^kit)) are periodic of period 
Tk. Because of the limited bandwidth of Sa^k^ the expan- 
sion in Fourier series of these two functions only involves 
frequencies 0, jr and — ^ of Sa^k- 



Our goal is to find filter g(z) producing a signal r(n) that 
coincides with a filtered version of one of the source sig- 
nals {sk)k=i....,K- This can be achieved by minimizing a cost 
function. In the following we investigate whether or not 
the Godard cost function is a good contrast function for 
mixtures containing BPSK signals. In a cyclostationary 
context and for a discrete time signal r, the Godard cost 
function is defined as 



J(r) =<E(|r(«)|2-l)' 



(4) 



In order to express J{r) in a more convenient way, we re- 
mark that r(n) can be written as 



r(n)^J]e^'''"'Hfk(z)]sk(n) 



(5) 



k=l 



where fkiz) is the transfer function fk(z) = g{ze ^"'^^'')hk(z). 
We denote by \\fk\\ the norm of filter fk(z) defined by 



J-l/2 



where Sf^^ie^'"^) represents the spectral density of signal 
(i/t(n))„ez- We finally define filter fk(z) and signal Sk(n) by 

(6) 



fkiz) = Skin) = [fk{z)]sk{n) 

WjkW 



If II/aII = 0, we put fkiz) - and Sk{n) = 0. It is clear that 
||_^|| = 1, and that < E|St(«)p >= 1. r(«) can be written as 



r(n)^Yj\\fk\\e^""'^'^Sk{n) 



(7) 



A-=l 



Note that when the excess bandwidth -yk is small, the cyclic 
correlation coefiicients at non-zero frequencies are clearly 
inferior to those corresponding to the zero cyclic frequency. 

We denote ak - k - I, . . . , K. Then, it is clear 

that the non zero cyclic frequencies of the discrete time 
signal are +ak] if moreover Sk is a BPSK signal, its non 
conjugate cyclic frequencies are 0, ±ak. From now on, we 
denote by / and Ic the set of all cyclic and non conjugate 
cyclic frequencies of y(«). We obtain immediately that 

. / = {0,(+a,)A.=i k] 

• 4 = {{26fk, 26 fk + ak)k=i....,K, Sk BPSK) 

In the following, we also denote /* the set of non zero cyclic 
frequencies of y(«). 

In order to extract one of the source signals, (y(n))„ez 
is filtered by a A^~inputs / 1-output filter g(z) to produce 
the 1-dimensional signal r{n) = [g(z)]y(«). It is straightfor- 
ward that this scalar signal r{n) has the same cyclic and 
non-conjugate cyclic frequencies as the received signal y(«). 



and coincides with a filtered version of one of the source 
signal (up to the term e^'™''/') if and only if the coeffi- 
cients {\\fk\\)k=\,....K satisfy \\fk\\ = 6{k - ko)\\fkJ\. We state the 
following result 

Proposition 1. The Godard cost function given by ^ 
can be expanded as 



J(r) = Y,mWk\f+J]Kh,,hMkf\mf-2j]\\fkf+i 

(8) 



k=l 



k,i^k. 



k=l 



where l(Sk^,Sk2) and /3{Sk) are defined respectively by 
2+Re 



and by 



(9) 



ael, 



'=-1,1 



/=- 1,0,1 
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Proof. We start by writing J{r) as 

J{r) =< E|r(n)|'* > -2 < E|r(«)|2 > +1 
Using the relation 

^rin)^ = C4(r(n)) + 2 (E|r(n)p)^ + |E(r-(n))f 

and the Parseval identities < (E|r(«)p)^ Y.aei l^r"'(0)F 

and < |E(r(n))-| >— ^aei l^c^r(0)Pj ^6 immediately get 
that 

7(r) C4(r(«)) > +2^ |/?;"^(0)P + ^ |/?[5^(0)p-27?J."'(0)+ 1 

ael ael, 

(11) 

Since (si:)ji=i,...,A: signals are independent we can write 

K 



< C4{r(n)) > ^ Yj "■^'"'^ ^ C4{sk{n)) > 

k=l 
K 



k=l 



where a represents one of the cyclic frequencies of r. For 
a -0, the last expression becomes Rf\Q) - 2f=i WfkW^ since 
we assumed that < E|Si(n)|^ >= Rf^\0) - 1. Furthermore, it 
is easily proved that if Oc is one of the non-conjugate cyclic 
frequencies of r, then the non-conjugate cyclic correlation 
coefficient of signal e^"'"^-^''Sk{n) at o-c frequency and at time 
lag coincides with R^"^^^^^''\Q). This implies that : 



fkfR^X''^'\0) 



k=l 



Using these various expressions in (jlip we obtain the an- 
nounced result. Notice that it is easy to establish that 
Pisk) is also given by 



(12) 



Note that, as shown in P(Sk) =< ElS^Cn)!"* >> 1. 

Expression ^ shows that J{r) is a function of both the 
norms {\\fk\\^)k=i,...,K, and the unit norm filters {fk{z))k=i,...K 
defined by Sk{n) - [fk(z)]sk(n), and that these 2 sets of 
parameters are independent. Minimizing J{r) with respect 
to g(z) is thus equivalent to minimizing ([8|) independently 
with respect to the norms (||/itlP)/t=i,...,A' and the unit norm 
filters (fk(z))k=\,...K- 

In the following we study the minimization of J{r) firstly 
when the different source signals do not have any non zero 
cyclic frequency in common nor any non-conjugate cyclic 
frequency in common and then we consider an opposite 
scenario where K BPSK signals share the same baud rate 
and the same carrier frequency. 



3. The source signals do not share the same cycHc 
and non conjugate cyclic frequencies 

We first study the behavior of J(r) when the source sig- 
nals do not share the same cyclic and non conjugate cyclic 
frequencies. This situation is likely to occur when the dif- 
ferent transmitters do not belong to the same network and 
it practically implies that / e {1 ...K} Ok + ai (i.e. 

Tk + Ti) and 6fk + 6fi {Afk + h.fi). In this context, the 
term lisk^,Sk2) reduces to the constant term 2, and J{r) is 
given by 



>ii'+i 



(13) 



k=l 



k=l 



We now study the conditions under which the minimum 
of J{r) is reached for a filter such that - 5(k - A;o)||/a„||. 
For this, we follow [l3| and we first fix the unit norm fil- 
ters (fk)k=i,...,K or equivalently the (JS(sk))k=i,....K coefiicients. 
Then, we consider the problem of minimizing J with re- 
spect only to the (||/A:lP)jt=i,...,A: • This is an easy task be- 
cause, as a function of the {\\fk\\^)k=\,...,K norms, J{r) has a 
simple expression which allows the following result to be 
derived 

Theorem 1. The minimum of J(r) w.r.t. (\\fk\\^)k=i,...,K is 
reached for sequences such that \\fk\\^ = ^(^ ~ ^o)II/aoII^ ^ 
certain ko index if and only if 

min B(sk) < 2 

k=\,....K 

and if this minimum is reached for the index ko. Moreover, 
the minimum value of J is equal to I — ^ ^ ^ . 

Corollary 1. // the sources do not share the same cyclic 
and non conjugate cyclic frequencies, the global minimiza- 
tion of the Godard cost function allows to extract all the 
source signals using a deflation approach if 

/3min,k^ min /3{sk) < 2, for each k ^ I,.. .,K (14) 

A,IIAII=i 

The proof of this theorem can be found in [l3|. It re- 
mains to check if condition holds. For circular lin- 
early modulated signals, (fn]) has been analytically proved 
in [lO|. In the case of BPSK signals, the following result 
can be proved using a similar approach. 

Proposition 2. Consider a BPSK signal with symbol pe- 
riod T and excess bandwidth < y < 1, and assume that 
the sampling period Tg does not belong to {T, De- 
note by K the kurtosis of the corresponding binary symbol 



sequence, k 
given by 



-2. Then, yS„ 



mm; 



inf 



ma) 



xPiUizmn)) IS 
(15) 



where <!)(/„) is defined by 



0(/„) = kT 



+ 2+4 



Tdt] 



'Tdt 



{k\fM?drf 



Moreover, if we define rjmin by ri„ 
then 



inf kT— 

/„<Er([-ii?,^]) ( |/„(f)l W 



< C4(S) >, 



(16) 



We give the proof of this result in the [Appendix A 



Remark 1. If e {T,^,^,^■}, the expression and, as 
a consequence, the value of jSmin, is different from il5\). 
Pmin, 0-s a function of Tg, is therefore a constant func- 
tion except in such points as {T,j,j,^} where it has a 
different value. We are therefore dealing with a discontin- 
uous function. In order to illustrate this point, we consider 
as an example the case where Tg = T . For this sampling 
rhythm to satisfy the condition of Shannon it is necessary 
and sufficient that the excess bandwidth factor be 0. 

Under these conditions, it is well known that Pmin equals 
1 and Tjmin — —2, while we will soon see that f3mm — 1-19 
and Tjjnin - —1-36 ifTe does not belong to {T, ^, j, =^]- For 
simplicity reasons, we prefer not to give the expressions of 
Pmin if Te ^ In cmy case, the probability of Tg 

being equal to one of these values is obviously null in a 
blind context. For this reason we suppose in the following 
that Te does not belong to {T,^,^,^]. 

As C if yi < yj, (USD 

implies that considered as a function of y, Pminij) is de- 
creasing. This observation allows us to make the following 
statement : 




(a) 



- N-PSK signals, 

- 64-QAM signals, k=-0.69 

- 256-QAM signals, tc=-0.6 
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(b) 



Proposition 3. Function j Pminij) is decreasing when 
J varies from io 1 . Consequently, jSmmiy) is strictly in- 
ferior to 2 for all y if and only if /3min(0) < 2. 

The main interest of proposition [3] is that if a function 
fait) 6 'F{[-^,4f]) (corresponding to y = 0), then the 
integrals 



Jffi Jr Jk 



vanish. This result is a direct application of the inequality 
of Parseval. The expression of ;8niin(0) is therefore 



;Smin(0) = min kT 



+ 2 + 



(17) 

It is easy to notice that jSmin does not depend of T and 
that the theoretical expressions (|TS|) and (fTOj) of y6,nin and 
'7min can be used in order to compute the numerical values 
of these functions for all the values of y € [0, 1] via the ap- 
proach proposed in [loj .Figure 1(a) gives a numerical rep- 
resentation of jSmin as a function of y in the case of BPSK 
signals. Moreover, we have found that T/min - 0.68/<'(l + y) 
and is equal to -1.36(1 + y) in the case of BPSK signals, 
since — —2. Figure 1(a) also confirms the decreasing 
nature of ySmin with respect to y, and the fact that for a 
BPSK modulated signal /3min < 2 for all y provided that 



Figure 1: ySmin as a function of y in the case of (a) BPSK signals and 
(b) cicular signals. 



does not belong to {T,T /2,T /3,2T /3}. If T, = T, as 
we have already mentioned, jSmin = 1; if equals one of 
the other possible values, we can directly verify that jSmin 
remains strictly inferior to 2. We can therefore enunciate 
the following result: 

Proposition 4. In the case of circular or BPSK trans- 
mitted signals, not sharing any non zero cyclic frequency 
nor any non conjugate cyclic frequency, the minimization 
of the constant modulus criterion, along with a deflation 
approach allows the extraction of all sources. 

Remark 2. Notice that the values of Pmin for a BPSK 
modulated signal are smaller than the ones we observe for 
linearly modulated circular signals which we represent in 



figure 1(b) . This means that if a BPSK modulated signal 



is mixed with circular modulated signals, the BPSK source 
will very often be the first one extracted when using a de- 
flation approach. 

Despite its undeniable importance, proposition [2] is not 
completely convincing as to the pertinence of the pro- 
posed approach. In practice, the search for filter g(z) — 
8(0z ' which extracts a source from the mixture is 
done by minimizing an estimator J{r) of J{r). Further- 
more, the minimization of 7(r) is carried out by means of 
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iterative algorithms such as the steepest descent or New- 
ton algorithms who are not guaranteed to converge toward 
the global minimum of J and may very well converge to- 
ward a local minimum instead. It is therefore necessary 
to verify that / does not have any non separating lo- 
cal minima. Under a technical assumption, the following 
result can be established 

Proposition 5. Assume that at least one of the functions 
fk — > /3([fk{z)]sic(n)) defined on the set of all unit norm fil- 
ters has no local minimum such that /3{[f^(z)]sic(n)) > 2. 
Then, the argument of each local minimum of the Godard 
cost function is a separating filter. 

Proof. We define the following quantities u = 



iZli WfkW^y^ and V, 
becomes 



ll/ill 



Expression (jl3p of J(r) then 



J(r) 



k=l h*k2 



2m' 



2' 

VA=1 



+ 1 



In the following we pose p - (fi{si), . . . ,li{sK)Y and we 
denote by T{\,^) the expression multiplying the term u^. It 
is clear that ^f^i = 1. Since Zk.^k. vlA.2 = (2^f=i ^'l^^ ' 
XfLi we obtain a simpler expression for T{y,P) 

K 

T{y,p)^2 + Y^vl(Ji{~su)-2) 



k=l 



J{r) is thus given by: 

J{r) = u'T{y,P) - 2ii^ + 1 

We consider a local minimum (f^iz), ■ . ., f^iz)V of J{r), and 
denote by «*, v,, /^*, the corresponding values of 

u,\,fk, Sk,p. It is easy to check that the point v, is a local 
minimum of the function v T{\,/3,). As at least one the 
coefficients (j6(sp - 2) is strictly negative, - S(k - ^o)v^^ 
where A:o is one of the index for which Pk^,* - 2 < (see 
e.g. i). This implies that = 6{k - ^o)ll/;jl, and 

that the local minimum . . . ,/a:,*(z) is a separating 

filter. It is difficult to check analytically whether or not it 
exists k for which fk — > P{[fk{zy\sk{n)) has no local minimum 
f^ such that P{[fk{z)]sk{n)) > 2. However, this condition 
probably holds because the steepest descent minimization 
algorithms of the functions ft — » P{[fk(z)]sk{n)) we have run 
always converge toward a point for which li([fk{z)]sk{n)) < 
2. 

In sum, the above results indicate that if the source 
signals do not share the same cyclic and non conjugate 
cyclic frequencies, then, the minimization of the Godard 
cost function allows to extract circular and BPSK source 
signals. In this context, it is therefore possible to separate 
the source signals without any knowledge of their cyclic 
and non conjugate cyclic frequencies. 



4. K BPSK sources sharing the same baud-rate 
and the same carrier frequency 

In this section, we consider the opposite situation, when 
all the source signals are BPSK signals with the same baud 
rate T, the same carrier frequency offset A/, and the same 
excess bandwidth y. We also denote by a and 6f the terms 
a - Te/T and 6f = AfT,,. Recall that the sampling rate 
is assumed not to belong to {T, T/2, T /3,2T /3}. 

4.1. Existence of spurious local minima for K — 2 and 
r = 

Our purpose is to support the conjecture that the Go- 
dard cost function has non separating local minima, and 
that the minimization algorithms often converge toward 
these spurious points. In order to justify this, we assume 
that the common excess bandwidth y of the 2 source sig- 
nals is equal to 0. In this context, the cyclic and non 
conjugate cyclic correlations coefficients at frequencies ±a 
are zero. Expression (|5]) of J{r) thus reduces to 



j{r) = Pimfit + m)\\f2t+ 



(18) 

211/1 ll'll/2ll' (2 + Re(7;^^_(0)/;[!']^(0)*)) - 2 (ii/iip + ii/jip) + 1 

where yS(i,) is given by 

PiSd =< c,{~sd > +2 + \R%{0)\' 

for / = 1,2. This expression is formally similar to the one 
of J in the case where the 2 sources are circular with a non 
zero excess bandwidth (see [l3l)i except that the cyclic 
correlation coefficients R^liQ) are replaced by 27?j"'(0). An 

analog of the condition |2i?j"''(0)| < 1, which plays an im- 
portant role in [loj . can also be proved true for the cyclic 
correlation coefficients R^^\{0), i.e. |/?fj,(0)| < 1. Consider- 
ing the definition of S,- in (|&]), we can write 



Jr 



As signal Sk is real valued, Sfl^, coincides with the spectrum 
Sfl^ of iji, and is an even function. Using the Schwartz 
inequality, we get immediately that ( 0)| < 1. It is 
therefore possible to use Theorem 2 of [lol established in 
the circular case to prove that if ySmin and /ymin defined in 
Proposition [5] verify 

f -3ySniin + 5 + ?7min > 

I 20S„iin - 1 )Ani„ - 4 ( 1 - i V206mi„ - 1) - 1 - 7?min)^ < 

(19) 

then, the argument of the global minimum of J{r) is a 
separating filter, and the minimum value of J{r) coincides 
with 1 - l/j8min- For r = 0, ySmin - 1-19, 7?,nin - -1.36, and 
it is easily checked that the 2 conditions above are satis- 
fied. The global minimization of J{r) therefore allows to 
separate the 2 BPSK signals. Moreover, 1 - l/y6min - 0.16. 
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However, J(r) may have non separating local minima, to- 
ward which a steepest descent minimization algorithm of 
J(r) often converges. In order to define these local minima, 
we denote by f^iz) one of the arguments of the global mini- 
mum of j8([/i(z)]ii(«)) over the set of unit norm filters with 
real coefficients. We denote by y6i,m,„ the corresponding 
minimum. It is easy to show that ySi,„„„ can be evaluated 
using Proposition [5J by minimizing the function over 
the real elements of 'F{[-l/2T, 1/27]) when y = 0. In 
these conditions, it can be shown that ySi^mm coincides with 
'7min + 3, i.e. that /3i,„„„ ^ 1.64. We now consider the unit 
norm filter with imaginary coefficients /jCz) = ifliz). It is 
clear that yS([^*(z)]i2(«)) coincides with P\,mi„. We finally 
define filters f*{z) for /= 1,2 by 



1 



(1 +y6i,„™)'/2 



Ynfiiz) 



(20) 



If r»(n) = [/i*(z)]«i(n) + [/2*fe)]*2(«), one can check that 
J(r,) - 1 - 2/(1 + Px^min) - 0.25. Although we have not 
been able to analytically prove these non separating points 
to be a local minimum of 7, we have observed that the 
steepest descent minimization algorithm of J{r) very often 
converges to one of these points rather than toward the ar- 
gument of the separating global minimum of J . To verify 



this, we present in Figure 2(a) an histogram of the values 
of J(r) at convergence of the steepest descent minimization 
algorithm. We used 1000 experiments, each corresponding 
to different randomly selected propagation channels, and 
we assumed the thermal noise to be negligible. The figure 
clearly shows that in more than half of the experiments 
the final value of J{r) corresponds to 1 - jT' ~ ^--^^ which 
is associated to a local minima rather than to the value of 
the global minimum of J which is 1 - = 1 - y^g - 0.16. 

In order to verify that the value 1 - does not corre- 
spond to a separating filter, we present in figure [2] an his- 
togram of the signal to interference and noise ratio (SINR) 
associated to the filters determined by minimizing J{r). We 
define the SINR as the ratio between the power of signal 
ri, representing the contribution of the extracted signal 
filtered by the extracting filter and the power of signal r2 
which represents the contribution of the other transmitted 
signal filtered by the same filter. It is clear that if the filter 
is perfectly adjusted then the SINR must equal +oo in the 
absence of thermal noise. The experiments we presented 
thus tend to confirm the fact that J{r) has non separat- 
ing local minima and that the steepest descent algorithm 
converges very often toward one of them. 

4-2. A new cost function 

A simple modification of the Godard cost function allows 
to overcome the aforementioned problems, provided that 
the most significant non-conjugate cyclic frequencies of the 
received signal are known or can be correctly estimated 
by the receiver. We recall that for a mixture of BPSK 
modulated signals sharing the same carrier frequency, the 
most significant cyclic frequency is 26 f . 



500 p 
450 - 
400 - 
, 350 - 
: 300 - 
' 250 - 
! 200 - 



0.12 0.14 



0.16 0.18 0.2 0.22 

CMA criterium values 



(a) 



10 15 20 25 30 35 

SINR 



(b) 



Figure 2: Histograms of (a) J{r) values and (b)SINR values, obtained 
after extracting a source from a mixture of 2 identical BPSK signal 
with 7 = 0. 



In the following, we assume that the carrier frequency 
offset 6f is known or correctly estimated at the receiver 
side, and consider the cost function J (r) defined by 



j\r)^J{r)-\R'^^/\0)\ 



(21) 



E(|r(n)p - l)" > - |< Mr^{n))e 



-2inn26f 



j'(r) is obtained by subtracting from J(r) the modulus 
square of the non conjugate cyclic correlation coefficient 
at time lag and at non conjugate cyclic frequency 26 f. 
Using the expression of we immediately obtain that 

K K 

J\r) = Y,P(~Sk)\\fkt+Yj /'(S*i>SaOII/aJ|'||/*,I|'-2 2 II/aII'+I 



i-=i 



k=l 



where the term I (Sj^, , is given by 



(22) 



2 + Re 



2 ^ /?^>(0)(/?<^;'(0)y+ ^ RU,(^)'^xjfiy 



'=-1,1 



/=-i,i 



(23) 
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and where (sk) is defined by 

< cdh) >+2 + 2Y, \Rt(0f + X l^cl'^O)!' (24) 



i=-i,i 



/=-i,i 



also equals 



(25) 



In order to give some insight on J , we first consider the 
case 7 = 0. Expression ((^ of /(r) therefore becomes 



J'(r) = Y^f^'ihWklf + 2 2 WfkfWfkf - 2^ ll/^ll' + 1 



k,*k2 



k=l 



Furthermore, now equals =< C4(S,) > +2. The 

expression of J (r) is thus similar to ((T^ , except that 
y6(S,) is now replaced by /3 (S,). It is easy to check that 
< C4(S,) >< 0, so that y6 (S,) < 2 for each /. Theorem [T] and 
Proposition[5]thus imply that the global minimum and the 
local minima of J are separating filters. This shows that 
the minimization of j'(r) allows to separate the K BPSK 
signals if y = 0. 

In order to extend this result to the more general case 
where y > 0, we now show that the argument of the mini- 
mum value of J'{r) corresponds to a separating filter. Con- 
trary to the case where the transmitted signals all had dif- 
ferent cyclic frequencies and different non conjugate cyclic 
frequencies, it is no longer possible to directly character- 
ize the global minimum of J (r) since its analytical form 
is too complex. Wc overcome this difficulty by using the 
following result stated and proved in [l^l : 

Proposition 6. Let m{r) be a positive function such that 
for any filtered version r(n) — [t(z)]s(n) we have 

J (r) > m(r) 

Assume that the infimum of m{r) is reached if and only 
if signal rin) coincides with a filtered version of one of 
the source signals. Let nin) - [fk(,,*(^)]sko(n) be one of the 
signals for which inff(j.) m(r) = m{rt). If m{rt) — J (r,), then 



infy (r) = y (r.) 



and the infimum is reached if and only if r{n) coincides 
with one of the r» specified above. 

In order to derive a function m{r) satisfying the condi- 
tions of Proposition [SJ we prove the following result. 

Proposition 7. The following inequality holds: 



Re 



(la) 



(0)* 



'=-1,1 



/=-l,l 



> -3/2 
(26) 



We give the proof of this result in [Appendix B Con- 
sider the function m(r) defined by 



m(r) = 



\k=\ 



\ki ^k2 



2 -2£||/,||2 + i 

(27) 



k=\ 



where we denote by yS . the quantity 



m!n,k 



with B . I - miHi 



\\fi_\\=i/3 (Sjt). Recall that the signals present 
in the analysed bandwidth are of the same nature and 
therefore aU 08^„^j^.)i^i,...,A: are equal. Since relation (pS)) is 
verified, it is clear that /'(Si, > 1/2. Moreover, the 
(J3' iSk))k=\,...,K are all greater than ySj^-,;,-,. This implies that 
for all r, j'(r) > m(r). We show that if < 1/2, then, 
the global minimum of m{r) is reached if all {\\fk\\)k=i,....K 
are null except for 1, i.e. if r(n) coincides with a filtered 
version of one of the sources. In order to establish this 
rcsuh, we pose = ^f^i \\fkf, = iif^, v = (vi, . . . , v^)^, 
and wc define f(v) as 

1 1 

^ A-=l ^ 

It is easy to verify that 

m(r) = u'^t(\) -2u^ + l 

and that the global minimum of m(r) is reached in a point 
(m,, V,) for which f(v,) is minimum and ul - The value 
of this minimum is then 1 — . To conclude it suffices to 
remark that if yS^^j^^ — 1/2 < 0, then the minimum of f(v) 
is reached if and only if all the components of v are null 
except for one who is equal to 1 , which corresponds to all 
WfkW being null except for one of them (0). Furthermore 
t(\t) is equal to ySj^-,;,-,, ul — and the minimum value 

of m{r) is 1 - In the following we denote by one 

of the index for which /3'^_^^^_^ = /3',„i„^k„j and by /a.„,, a unit 
norm fihcr for which = {[fkoAz)]skoin)), and we 

pose fko,*{z) - Utfko,*(z). The minimum of m(r) is reached 
if r,(n) - [/io,,(z)]i/io(")j and J (r,) coincides with m(r,) - 
l-^p—. Proposition|6]then states that the global minimum 

of J is reached only if r(n) is a filtered version of Skg(n). 
We have thus established the following result 

Proposition 8. If P^^^ ^ ^/2, then the minimization of 
J (r) allows the extraction of one of the sources from the 
mixture. 

Wc must now verify whether the condition B . < 1/2 

^ " mill ' 

is satisfied or not. Following the same reasoning as in the 
case of j6min, we can easily adapt proposition [2] by simply 
replacing the expression ((T5|) with 



inf 



(28) 



where <1) (fa) is defined as 

<1> (fa) ^kT^ — +2 + 4 



(29) 



The expression of ^'{fa) is obtained directly by subtract- 
ing from the expression of <!'(/„) ([TB)) the term due to the 
square modulus of the non conjugate cyclic coefficient of 
s{n) at the non conjugate cyclic frequency 0. As in the 
case of ySmin, this result implies that fi . is a decreasing 
function of the excess bandwidth factor y. We can thus 
formulate the following statement: 

Proposition 9. The function j — > y6,„j„(T) is decreasing 
when y varies from Q tol. Consequently, j8„j„(7) is strictly 
inferior to 1/2 for all values ofy if and only ifPminiO) < 



The expression of 
the one of j6min(0) ^ 



(0) can be deduced directly from 



^mi„(0) 



mm 



kT 



+ 2 = 77n,in+2 (30) 



with 77min given by equation ((T5)) . 

Recall that we can numerically evaluate the values of 
^min Imin for all excess bandwidth factor y e [0, 1]. 
Particularly for an excess bandwidth factor of 0, ?7min — 
-1.36 and . = 0.64 > 1/2. In order to verify the ex- 

" min ' 

istence of some 0^^^ values smaller than 1/2 we present 
in figure [3] the graph of 0^^^{y) for all excess bandwidth 
factor ye [0, 1]. The figure confirms the decreasing nature 
of y6,-,-,ij-, with respect to y and shows that yS^j^ < 1 /2 as 
soon as y > 0.1. Consequently, we are sure to separate 
the BPSK sources using the minimization of j'(r) if their 
common excess bandwidth factor is superior to 0.1. 




Figure 3: 



as a function of y for a BPSK signal. 



When the excess bandwidth factor is inferior to 0.1, the 
inequality j'(r) > m{r) does not allow any conclusion to be 
drawn as to the global minimum of J (r). However, in such 



cases, we can consider the approach used in [loj in the case 
of circular signals and inequality (^5)) . After some algebra, 
we can prove that if j0 . > 1/2, then the global minimum 
of J (r) is reached for filters which allow the extraction of 
one of the sources, if the following 2 sufficient conditions 
are met. 



77„,in + 3 - + 1)08„ 



i)>0 



(K - 1) (2 - ^_^{K(0,,,^ - i) - (77„,i„ + I)) <0 

^(31) 

We give the proof of these conditions in [Appendix C| We 
can easily verify that these conditions hold for ye [0,0.1] 
if the number of sources K is inferior to 10, which is very 
satisfying in the considered context. 



5. The case of general mixtures 

5.1. Generalisation of J (r) 

The results obtained in the case of a mixture of BPSK 
signals sharing the same characteristics can be extended 
to more general mixtures of circular linearly modulated 
signals and BPSK signals. The logic behind the definition 
of j'{r) is to subtract from J{r) the square modulus of the 
non conjugate cyclic correlation coefficients at time lag 
and at the non conjugate cyclic frequencies {26fk, Sk BPSK}. 
These frequencies are called in the following the significant 
non conjugate cyclic frequencies of the received signal, and 
we denote by Ic^s this set. The definition of j'{r) thus 
becomes 

/(r) = 7(r)- 2 lijf'^'W (32) 
=< E(|r(«)|2 -^f >- J]\< E(r^(n))e-'''""'' >\' 



We assume that the mixture contains L groups of (_Ki)i^\^_^L 
BPSK signals sharing the same characteristics (symbol pe- 
riod, carrier frequency, excess bandwith) and linearly mod- 
ulated circular source signals whose symbol period differ 
from those of the BPSK signals. If source k is circular, 
then it holds that ■ , - Bmi„k > I > . Therefore, it 
is easy to check that J (r) > m{r) where m{r) is still de- 
fined by ([27]) . Proposition [8] thus implies that if the excess 
bandwith of the BPSK signals are greater than 0.1, then 
the minimization of J allows to extract all the BPSK sig- 
nals. The case where some of these excess bandwiths are 
less than 0.1 is more difficult, but could be addressed using 
the previous approach. We just mention that if the cyclic 
and non conjugate cyclic frequencies of the sources are 
pairwise different, then the minimization of J still allows 
to extract the K sources whatever their excess bandwiths. 
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In effect, J (r) can be expressed as 

K K 

Hr) - 2y8'(*^)||/*ir + 2 ^ \\fk,f\\fk,f-2Y,\\fk\\- + 1 

k=l hi^ki k=V 

(33) 

If the source k is circular {sk) coincides with fi(sk) (fTU)) 
while for a BPSK source is defined by The 

expression of J (r) is therefore similar to that of J{r), and 
thus all results obtained in section |3] remain valid since for 
all Pkmin — Pkjnm < 2. The modification of J proposed in 
order to solve the problems generated by mixtures of non 
circular sources of the same nature thus does not modify 
the results obtained in the context of circular or non 
circular sources having different cyclic and non conjugate 
cyclic frequencies. 

5.2. Frequency offset estimation 

The use of J requires of course the correct estimation 
of the significant non conjugate cyclic frequencies of the 
received signal prior the source separation. Fortunately, 
this is a much easier task than the estimation of the baud 
rates, because the non conjugate cyclic correlation coef- 
ficients of the received signal at twice the frequency off- 
sets are not affected by possible low excess bandwidths 
of the source signals. A simple detection technique based 
on the examination of the modulus of the periodogram of 
the signal {ym{n + T)y,„(«))«EZ (see for example 0) may be 
successfully used. We also notice that if the estimation 
algorithm detects not only the significative non-conjugate 
cyclic frequencies {25/^., BPSK}, but some non significa- 
tive conjugate cyclic frequencies such as 26 fk^ + au^ or 
25fko ~^koj then the behaviour of function J is even better 
because /3'{skJ defined in principle by (^5)) is replaced by 

mj-\Rn/-\i^^\' ormj-\i^/-\ii::Y is 

thus lower than what is predicated by Figure [3l The suf- 
ficient condition , > I is thus less restrictive than in 
the case where 25/^;^ + ai^^ and 25/^^ — ai^^ are not detected. 

6. Simulations 

6.1. Implementation of the deflation approach 

In order to introduce the deflation approach we have 
implemented, we consider g' the extracting filter obtained 
by minimizing the cost function {J{r) or J (r)). We denote 

by 

him) = {g\z)]y(m) 

an estimator of a filtered version of one of the source sig- 
nals. The deflation approach consists in subtracting the 
contribution of this particular source from the observed 
signal y(m). As a result a new signal y'^*(m) is formed con- 
taining only the contributions of the other sources. We 
can then run the extraction algorithm on y^^\m) in order 
to determine a new filter g^^\z) for which, signal 

him) = [g<'^(z)]y^2)(,„) 



represents an estimator of a filtered version of a second 
source. In practice, the first subtraction is not perfect 
and ri(m) may contain residual filtered versions of the re- 
maining sources. This can render the convolutive mixture 
defined by y^im) more difficult to inverse than the origi- 
nal one. It is therefore reasonable to try to go back to 
the originally mixture yitn), and apply the extraction al- 
gorithm initialised with a filter close enough to g*^'(z), the 
filter that allows the extraction of the second source. This 
initial filter, denoted by gfliz) = Z^^q S^tj^^'^ '^^ obtained 
by minimizing with respect to g|^?, = istifl' SMt,D ^^e 
quadratic criterion 

J M-l 

m-O 

This initialization, proposed in jl4| . allows the extraction 
of the second source from the original mixture y to be 
achieved with better performance. 

6.2. Simulations parameters 

The experimental results we present in the following 
were obtained in the context of blind separation of a con- 
volutive mixture of K - 3 equal power BPSK modulated 
signals, observed by a receiver equipped with a circular ar- 
ray of N - 5 sensors distanced from one another by half a 
wavelength. All sources have the same excess bandwidth 
factor y = 0.5. 

The propagation channels are multi path and affected 
by a Rayleigh fading. An arbitrary path (k) is charac- 
terized by its delay Tk, elevation ^i,, azimuth 6k and at- 
tenuation Af;. We consider the ETSI channels BUx, TUx, 
HTx, RAx. For each experiment, the arrival angles on 
the different paths {(p and 9) of the signals are randomly 
chosen inside [-;r/2, ;r/2] and [-n,7t] respectively. The dif- 
ferent complex amplitudes on each path are also randomly 
chosen for each experiment. Generating different channel 
characteristics from one experiment to another enables us 
to have statistically significant results. We suppose that 
the central frequency of the receive filter of the receiver is 
fo - IGHz and that the received signal is corrupted by a 
white, additive complex gaussian noise with power spec- 
tral density A^o- The signal to noise ratio per source signal 
^ is equal to 20 dB. We have considered two opposite 
scenarios : 

• all BPSK signal have the same symbol period T - 
3.6fis and same frequency offsets {6f) 

• the BPSK signal have different symbol periods {Ti = 
3.4fj.s, T2 - 3.6jj^s, T2 - 3.9fis ) and different frequency 
offsets {6fi +5f2 +5f^) 

In both cases, the sampling period is equal to ^g, and 
the carrier frequency offsets are randomly chosen on each 
trial such that the generated signals satisfy the sampling 
theorem. We also considered different observation dura- 
tions r„fe, = 2000r, Tob, = lOOOr and r„fe, = 500r for the 
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initial received signal used to compute the separating fil- 
ters and a longer observation duration of Tperf - 20000r 
for the performance analysis. For each possible type of 
mixture we considered 1000 independent experiments. 

6.3. Numerical results 

We chose two metrics of performance for our separating 
method : the signal to interference plus noise ratio (SINR) 
at the output of the separating filter g and the symbol 
error rate (SER) computed after applying a blind CMA 
fractional equalizer, supposed to know the baud rates and 
the carrier frequency offsets of the sources, to the sepa- 
rated signal. In order to compare the different separating 
algorithms we consider the number of experiments where 
we obtain a SER inferior to 10"^. 

Moreover, since the channels are randomly selected from 
one experiment to another, we need a reference measure 
of the difficulty of the separation problem. We chose to 
compute, for each source k the performances obtained 
in a non blind context with the minimum mean square 
estimator (MMSE). The Wiener filter g^^„^,(z) obtained 
with this method is a finite impulse response filter of the 
same size as g. This filter is chosen non causal, and its 
coefficients are estimated from the samples of the received 

signal (y('«))m=o m-i and those of the transmitted signal 

(S;t(m))m=o,...,M-i as if the receiver worked with a learning 
sequence of M samples. The performances provided by 
this filter thus represent an upper bound as to what we 
could achieve in a blind context. 



Table 1(a) contains the results associated with the first 



scenario. Notice that the number of times where the 
SER corresponding to the separation method based on the 
CMA algorithm is inferior to 10"^ is smaller than the one 
corresponding to the separation with the modified CMA 
criterion. This is due to the large number of cases where 
the CMA algorithm does not correctly extract the sources 
from the mixture. Contrariwise, the modified CMA algo- 
rithm succeeds in extracting one source from the mixture. 
This phenomenon is visible in figure U] where we present 
the histograms of the SINR obtained after the extraction 
of one source from the mixture using the CMA, modified 
CMA and MMSE methods, when considering BUx type 
communication channels and a duration of observation of 
2000r. It is easy to see that in an important number of 
cases the SINR values corresponding to the CMA method 
are close to dB meaning that no source was correctly 
extracted. The modified CMA algorithm significantly re- 
duces the number of unsuccessful extractions and its per- 
formance is close to that of the Wiener filter (MMSE). This 
phenomenon can also be observed on the results obtained 
on the other channels and when the duration of observa- 
tion is smaller. 

The results obtained in the second scenario are pre- 
sented in table 1(b) Notice that in this case the per- 
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CMA 
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Figure 4: Histogram of SINR obtained after extracting one source 
from a mixture of 3 identical BPSK signals, sent over a BUx channal 
and observed for a duration Tuts = 20007'. 



that of the modified CMA method but does not gener- 
ally surpass it. This confirms the good behaviour of the 
CMA algorithm when separating mixtures of signals with 
different characteristics all the while showing that the use 
of the modified CMA algorithm in such cases could bring 
some improvement. Generally, the performance of the two 
blind separation methods are close to that of the MMSE 
method. 



formance of the CMA extraction method is very close to 



7. Conclusion 

We investigated the separation of convolutive mixtures 
of second order circular linearly modulated signals and 
BPSK signals in the context of passive listening. We con- 
sidered only deflation approaches coupled with the min- 
imization of the CMA cost function. We proved that if 
the different source signals do not share the same cyclic 
and non conjugate cyclic frequencies, the minimization of 
the CMA cost function ensures the extraction of a filtered 
version of one of the source signals. We have also shown 
that in this case and under a condition which is always 
verified in practice, all the local minima of the CMA cri- 
terion are separating points. This result is no longer true 
when mixtures of BPSK signals sharing the same baud 
rate and carrier frequency are considered. In this case we 
have shown the existence of non separating local minima 
of the CMA cost function that prove to be quite attrac- 
tive. A modification of the CMA criterion was proposed, 
based on the knowledge of the most significant non con- 
jugate cyclic frequencies of the received signal. Moreover, 
the minimization of this new criterion was also proved to 
be a reliable approach in a much more general context. 
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Table 1: Percentage of TES < 10 ^ for a mixture of 3 BPSK signals whose cyclic and non conjugate cyclic frequencies are (a) identical and 
(b) all different. 



(a) 



No.symboles 


2000 


1000 


500 


BUx: 


BPSK 


BPSK 


BPSK 


BPSK 


BPSK 


BPSK 


BPSK 


BPSK 


BPSK 


CMA 


84% 


83% 


83.8% 


82.6% 


83.2% 


81.8% 


83.2%. 


81.8% 


84.7% 


CMAm 


99.9% 


100% 


100% 


100% 


100% 


99.8% 


97.4% 


95% 


97.2% 


MMSE 


100% 


99.8% 


100% 


99.9% 


100% 


99.9% 


99.8% 


100% 


100% 


TUx: 


BPSK 


BPSK 


BPSK 


BPSK 


BPSK 


BPSK 


BPSK 


BPSK 


BPSK 


CMA 


88.9% 


87% 


89.2% 


89.4% 


85.2% 


86.1% 


86.8% 


86.8% 


87.6% 


CMAm 


99.9% 


100% 


100% 


99.8% 


99.8% 


99.6% 


95.7% 


95.3% 


94.5% 


MMSE 


100% 


99.8% 


100% 


100% 


100% 


99.9% 


100% 


100% 


100% 


HTx; 


BPSK 


BPSK 


BPSK 


BPSK 


BPSK 


BPSK 


BPSK 


BPSK 


BPSK 


CMA 


89% 


87.2% 


86% 


87.6% 


88% 


88.4% 


87.5% 


87.7% 


86.1% 


CMAm 


99.7% 


99.8% 


100% 


99.2% 


99.5% 


99.4% 


91.5% 


91.5% 


91.4% 


MMSE 


99.9%. 


100% 


100% 


100% 


100% 


100% 


100% 


100% 


100% 


RAx: 


BPSK 


BPSK 


BPSK 


BPSK 


BPSK 


BPSK 


BPSK 


BPSK 


BPSK 


CMA 


78% 


79.7% 


79% 


78.9% 


81.3% 


79.6% 


81% 


81.1% 


80.3% 


CMAm 


100% 


99.9% 


99.9% 


99.4% 


99.1% 


99.2% 


93.6% 


94.9% 


92.9% 


MMSE 


100% 


99.9% 


99.9% 


100% 


100% 


100% 


100% 


100% 


99.9% 


(b) 


No.symboles 


2000 


1000 


500 


BUx: 


BPSK 


BPSK 


BPSK 


BPSK 


BPSK 


BPSK 


BPSK 


BPSK 


BPSK 


CMA 


100% 


99.9% 


99.8% 


99.8% 


99.7% 


99.8% 


99.5%. 


99.4% 


99.6% 


CMAm 


100% 


99.9% 


100% 


100% 


99.9% 


99.9% 


99.7% 


99.5% 


99.8% 


MMSE 


100% 


100% 


100% 


100% 


100% 


100% 


100% 


100% 


100%. 


TUx: 


BPSK 


BPSK 


BPSK 


BPSK 


BPSK 


BPSK 


BPSK 


BPSK 


BPSK 


CMA 


99.6% 


99.5% 


99.7% 


99.3% 


99.5% 


99.6% 


99.7% 


99.4% 


99.2% 


CMAm 


100% 


99.9% 


99.7% 


99.7% 


99.7% 


99.6% 


98.1% 


98.9% 


98.4% 


MMSE 


100% 


99.9% 


100% 


100% 


100% 


100% 


100% 


100% 


100% 


HTx: 


BPSK 


BPSK 


BPSK 


BPSK 


BPSK 


BPSK 


BPSK 


BPSK 


BPSK 


CMA 


98.8% 


98.7% 


98.8% 


98% 


97.8% 


98.1% 


96.7% 


96% 


96% 


CMAm 


100% 


99.8% 


100% 


99.5% 


99.4% 


99.4% 


98% 


98.2% 


97.7% 


MMSE 


100% 


100% 


100% 


100% 


100% 


100% 


100% 


100% 


100% 


RAx: 


BPSK 


BPSK 


BPSK 


BPSK 


BPSK 


BPSK 


BPSK 


BPSK 


BPSK 


CMA 


98.3%. 


98.9% 


98.3% 


99.2% 


98.9% 


98.5% 


99% 


98.8%. 


98.7% 


CMAm 


99.8% 


99.8% 


100% 


98.9% 


99% 


99% 


98.4% 


98.2% 


98.3% 


MMSE 


99.9% 


100% 


100% 


100% 


100% 


100% 


100% 


100% 


100% 



Appendix A. Proof of Proposition [2] 

Although the proof of this proposition is very similar 
to the one in [loj . we provide it in order to make the pa- 
per reasonably self-contained. For simplicity reasons, we 
assume the carrier frequency offsets to be 0. This assump- 
tion docs not reduces the generality of the results. 

Proposition 10. Suppose that is not a multiple ofT/2 
(this automatically holds from the hypothesis) and that 
Te < T/{1 + y) (this holds since^ holds). Let f{z) he a 
transfer function for which \\f\\ < oo. If g„{t) denotes the 
shaping filter of signal Sa(t), then the function fa(v) defined 
as 

/«(v) = /(e^^^^Ot^M, Vy e R (A.l) 

vanishes outside S — [— ^] and belongs to the space 
'F(JB). Let fa{t) he its inverse Fourier transform in the L^- 
sense. For every, t we define the continuous-time signal 



ira(t))tes. as 

rait) ^Y^a{j)Ut - jT) 

Then, the discrete-time signal r{n) — [f{z)]s{n) coincides 
with the discrete-time signal ra(nTe). 

This result is proved in [llj when the filter f{z) has 
a summable impulse response and in [loj when ||/|| < 
CO. Using proposition [TUl we have V/(z), 3fait) such 
that < E\r(n)f >=< E\ra(nT,)f > and < E{r(n)2} >=< 

E|ra(«re)^| >. Considering the time average of the 

Fourier series expansion of E|r„(f)|^ and E{r„(f)^}, when 

{r,f,f,f } wcget 

< E |r(n)|2 >=< E \r„{nT,)\^ > 

^RfJiO)^^ f IfMfdt (A.2) 



sion (l26l) as 



E{r(nf) >=<E{r,(«r,)'} > 



?(0) 



R'ciiO)-^^ r fa 



(tfdt 



(A.3) 



Furthermore 



¥.\ra{nT,fe^^'"Tdt. 



^ Jo 

■'Jo Jr 



Jk 



(0) 



/?^/(0) 



/l^/«(v)/«(y-i)*flfv 

2T 

fty fa(y)fa(^ - y)dv 

IT 

fi \fa{y)Uv 

IT 

IT 

fi\fa{y)Uv 



(B.l) 



(tfe^^'''Tdt 



In order to simplify the notations, we denote by and 
(A. 4) respectively B:^ the intervals B^ - ["^,"57"] and B^ — 
1^—^, -^]- It is straightforward that 



itr 

J " |/„(v)|' <iv > |/„(v)|' fl'v + j^^ \Uv)f dv (B.2) 



A similar reasoning can be carried on for the 4-th order 
cumulant of signal rg. function Cii{ra{t)) can be written as 
CA{ra{t)) - kYjii ~ w^l"* where k is the 4-th order cumu- 
lant of the transmitted symbol sequence. This function is 

periodic of period T and, due to the limited bandwidth of at ii ^ ii r i- r / \r / i\« 1 

7,, r , „ -TT^ - • Motice that the functions v — > ja(v)Ja{y - 7) and v 

hlter fa, has at most 7 cyclic frequencies. Its Fourier series 

expansion therefore is : 



Because of the conditions imposed on T^, the terms 
< e^'"''^ > arc zero if ^ 0, which means that: 



(A.5) 



faiy)faij — y) are zero unless v e B^ and that function 
V — » fa{y)fa{—y — 7) is also zcro outside of By. This means 
that 

Hz 

f" /«(v)/«(y-^)Vv= r /„(v)/«(v-^)Vy 

J- it? r r 

1+2 

J- it? r J fit r 



< C4(r(n)) CO = - ( \fa{t)\ dt 



hi 

My)M-y-l:)dy^ f iMfa 



(-y--)dv (B.3) 



where cq is the constant value of the Fourier series expan- Using the inequality of Schwartz we immediately obtain 
sion of < C4{ra(t)) >. 

Applying (TO]) . ((X3)) . (|X4)) and (|X5|) to signal 
s = s{n) and recalling that 



ii/ir 



/-1/2 

|/(e2,.v)|-5f«'v=<E|r/(z)l^(«)|2> 

J-1/2 



E|r/(z)li.(«re)|^ >= - lUtTdt (A.6) 



we evaluate the terms in expression (|15l) and find the ex- 
pression of ^ifa)- 



Appendix B. Proof of Proposition [7] 



In the following, we consider once more filter /a(v) de- 
fined in (|A.ip We begin by expressing the cyclic and non 
conjugate cyclic correlation coefficients involved in expres- 



r lie V ( r \ 1 P V 

j^^ h(v)!a(v- jTdv < \ j^^ |/;(v)f rfv ; j^^ |/„(y- -)| dv\ 

( r W r V'^ 

= 1^ J^^ |/;(v)|' ^/vj I |/<,(v)f dvj (B.4) 
j^^ /«(v)/„(^ - v)rfv < I j^^ |/„(v)|' rfvj j^^ 1 - v)|' rfy 

= r |/.(v)|'rfv 



(B.5) 



(B.6) 
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3^), we get 



Using (jB.ip , (jB.2|) , (jB.Sp and the means inequality {ab < 



\Rf\0)\ < 



(B.7) 



Since ^^""^(0) = (Rf\0)y we also get |^^"\0)| < i. 
We define the following quantities 

|2 



|/oi-,(v)| dv,yi= \fl,Xy)\-dy,ie {1,2} 



BZ 



= 2 (<;(0)4j(0)- + Rf;\o)Rf^(or) 



(0)/?r (0)* +7?pxo)/?r'(o)* 

where ki,k2 & {1,2} , ki ^ k2. 

With this notations and using (|B.ip . (|B.2|) and (|B.5p for 

signals Sk^ et 5^., , we get: 



3(a) 



< 4— + 



Xl 



X2 



x\+y\X2+y2 xi+yiX2+y2 x\+yiX2+y2 



(xi + y\){x2 + ya) (-^i + y\){x2 + yi) 
Using Schwartz inequality 



(B. 



^IX\X2 



^yiyi 



= y/(xi +yi)(x2+y2) (B.9) 



and since -j-^ < 5, we get that|f| < 4- 

Appendix C. Proof of conditions (|3ip 

Wc pose A(s,) = ^|4\Rf^\Q)f + (0)f + \Rt'\0)f . 

From (|^ . wc easily get A(sk) < 

Considering the expression of StJ it is easy to 
prove (using Schwartz inequality) that 

l{h„h,)>2-AiSkM(Sh) (C.l) 

Function j'{r) given by ([^ is therefore lower bounded 
by: 



fir) > £||/;i||";6'(S^)+ 2 



All! WfkiW 



A better lower bound for J (r) can be found by choosing 
a better lower bound of / (sk,, 5^) ^^'^ thus a good upper 
bound for A(sii^)A{sic2)- We first state the following obvious 
result. 

Lemma 1. Let Sa,k(t) - TjI '^i,kga,k(t - ITk) be one of the 

source signals and let Sk{n) be the discrete time signal ob- 
tained by sampling Sa^kit) at a rate ofT^. Consider a unit 
norm filter fk{z) and signal Sk{n) - \^fk{z}\sk{n). We con- 
sider an element /I, 6 (0; -Jl). If A{sk) > /l*, then: 



with 



where, /3 (sk) and 77 
tively. 



/?! = T]^^n + 2 + AI (C.2) 

mm are given by il5\) and il6]) respec- 



In the following we denote by "TCi the set of index k for 
which A{sk) < Af and by 'K2 the set of index k for which 
A{sk) > At. It is easy to see that this two sets are disjoint 
and that there union is the set of all source index TC = 
{I...K}. 



Observe that if ki,k2 e "TCi, then from (jC.lj) / (Sa, , %) > 
2-AI. If ki,k2 G 'K2, then s^.^) > and if ki e 'Ki and 

k2 € 'K2, then /'(S^i, Sa-,) > 2 - ^/l*. 

Under this observations we can give the expression of a 
better lower bound of J (r) : 

'«('•) =A„i„ 2 ii-^'^ii' + (2 " Z ii-^*' "'11^*^11' 



kii^kiMM^'Ki 



+ 2(2- J -At) ^ II/aJI'II/a, 



kie'K,,k2e'K2 



■2 2 ||/aI|'-2 2 11/^11' + 1 



(C.3) 



We consider a vector (||/a||)a=i,...,a' at which the gradient of 
m(r) is zero. Then, it is easily shown that 

• For each k e'Ki, ||/a|P is cither equal to or to a fixed 
value (independent of k) denoted by fi . We denote by 
Pi the number of non-zero terms. 

• For each k e'K2, \\fk\\~ is either equal to or to a fixed 
value (independent of k) denoted by fa. We denote by 
P2 the number of non-zero terms. 



and it can 



A=l 



ki*k2 



If f 2 = then fi is given by fi - --^ 

{2-A{hM(^k.)) . A„i„+(2-«;)m-i) 

be shown that the corresponding value of m(r) increases 



-2 2 II/* 



when Pi increases if 



k=l 



A„in^(2-4) 



(C.4) 
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If Pi =0 then t2 is given by fi - -y— — 7- and it can be 
shown that the corresponding value of m(r) decreases when 

P2 increases. 

If f 1 > 1 and P2> I, then fi and t2 are given by 



(2- yl^_A,)P2~(fi, + ^(P2-l)) 



(2 - ^jA.fPiP. - (fimin + (2 - Al){Pi - 1))08. + kP2 - D) 

(2 - ^^.)Pl - CSniin + (2 - ij)!/-! - 1)) 
(2 - ^jA.fPlP, - 08n,in + (2 - iJ)(Pi - 1))08. + {(P2 - D) 



('?min + 2; 77min + j). This mcans that wc should chose 
the smallest value for yS, which also verifies the first of 
conditions (|C.10p . 

For an excess bandwidth ye [0,0.1] a suitable value is 
j8, - Kifiram - \) + \- Replacing this value in equations 
(jC.lOp . we obtain conditions ([31]) . 
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(C.5) 

The corresponding value of m{r) is 

2(2 - ^A.)PiP2 - (Plifi. + i(P2 - D) + P2Ce,nin + (2 " AI){Pi ~ 1))) 

m = 1 p 
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for(Pi,P2)^(l,0) 



(C.7) 



In the following, we derive sufficient conditions for which 
((UJ)) holds. Wc first consider the case (^1,^2) = iO,K), 
and obtain the following condition 



/3, + \{K-l) 
K 



When (Pi,P2) = (1,^^- 1) condition ((UJ)) becomes 
(2- ^\)\K-\) 



>0 ■ 

— ^^min 



(C.8) 
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Replacing A, from jUJ) in jUH), dSll) and we 
obtain 

I Knin(fi' + j(K - 2)) -iK-l){2- ^|06. - 2 - 77„,in)) < 
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The minimum of J (r) is reached if all but one ||/^.|| norms 
are if conditions IC.IOI are simultaneously verified. This 

happens for a well-chosen value of /I* e (0, -^) and im- 

pHcitly of /?* € (?7i„in + 2; ?7min + 5). 

Remark that, as a function of j6,, 



^2 - -y/yS, - 2 - ?7min j is decreasing when yS* 
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